IRIX Base Documentation 2001 May/SGI IRIX Base Documentation 2001 May.iso/usr/share/catman/p_man/cat3/complib/BANDV.z/BANDV

Text File | 1998-10-30 | 6.8 KB | 133 lines

____BBBBAAAANNNNDDDDVVVV((((3333FFFF))))                                                          ____BBBBAAAANNNNDDDDVVVV((((3333FFFF))))



NNNNAAAAMMMMEEEE
     BANDV, SBANDV   -  EISPACK routine.  This subroutine finds those
     eigenvectors of a REAL SYMMETRIC BAND matrix corresponding to specified
     eigenvalues, using inverse iteration.  The subroutine may also be used to
     solve systems of linear equations with a symmetric or non-symmetric band
     coefficient matrix.


SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
          ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  bbbbaaaannnnddddvvvv((((nnnnmmmm,,,, nnnn,,,, mmmmbbbbwwww,,,, aaaa,,,, eeee22221111,,,, mmmm,,,, wwww,,,, zzzz,,,, iiiieeeerrrrrrrr,,,, nnnnvvvv,,,, rrrrvvvv,,,, rrrrvvvv6666))))
          iiiinnnntttteeeeggggeeeerrrr         nnnnmmmm,,,, nnnn,,,, mmmmbbbbwwww,,,, mmmm,,,, nnnnvvvv,,,, iiiieeeerrrrrrrr
          ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn eeee22221111
          ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn aaaa((((nnnnmmmm,,,,mmmmbbbbwwww)))),,,, wwww((((mmmm)))),,,, zzzz((((nnnnmmmm,,,,mmmm)))),,,, rrrrvvvv((((nnnnvvvv)))),,,, rrrrvvvv6666((((nnnn))))

          ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssbbbbaaaannnnddddvvvv((((nnnnmmmm,,,, nnnn,,,, mmmmbbbbwwww,,,, aaaa,,,, eeee22221111,,,, mmmm,,,, wwww,,,, zzzz,,,, iiiieeeerrrrrrrr,,,, nnnnvvvv,,,, rrrrvvvv,,,, rrrrvvvv6666))))
          iiiinnnntttteeeeggggeeeerrrr         nnnnmmmm,,,, nnnn,,,, mmmmbbbbwwww,,,, mmmm,,,, nnnnvvvv,,,, iiiieeeerrrrrrrr
          rrrreeeeaaaallll             eeee22221111
          rrrreeeeaaaallll             aaaa((((nnnnmmmm,,,,mmmmbbbbwwww)))),,,, wwww((((mmmm)))),,,, zzzz((((nnnnmmmm,,,,mmmm)))),,,, rrrrvvvv((((nnnnvvvv)))),,,, rrrrvvvv6666((((nnnn))))



DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
     On INPUT
     NNNNMMMM must be set to the row dimension of two-dimensional array parameters
     as declared in the calling program dimension statement.

     NNNN is the order of the matrix.

     MMMMBBBBWWWW is the number of columns of the array A used to store the band
     matrix.  If the matrix is symmetric, MBW is its (half) band width,
     denoted MB and defined as the number of adjacent diagonals, including the
     principal diagonal, required to specify the non-zero portion of the lower
     triangle of the matrix.  If the subroutine is being used to solve systems
     of linear equations and the coefficient matrix is not symmetric, it must
     however have the same number of adjacent diagonals above the main
     diagonal as below, and in this case, MBW=2*MB-1.

     AAAA contains the lower triangle of the symmetric band input matrix stored
     as an N by MB array.  Its lowest subdiagonal is stored in the last N+1-MB
     positions of the first column, its next subdiagonal in the last N+2-MB
     positions of the second column, further subdiagonals similarly, and
     finally its principal diagonal in the N positions of column MB.  If the
     subroutine is being used to solve systems of linear equations and the
     coefficient matrix is not symmetric, A is N by 2*MB-1 instead with lower
     triangle as above and with its first superdiagonal stored in the first
     N-1 positions of column MB+1, its second superdiagonal in the first N-2
     positions of column MB+2, further superdiagonals similarly, and finally
     its highest superdiagonal in the first N+1-MB positions of the last
     column.  Contents of storages not part of the matrix are arbitrary.

     EEEE22221111 specifies the ordering of the eigenvalues and contains
        0.0E0 if the eigenvalues are in ascending order, or



                                                                        PPPPaaaaggggeeee 1111






____BBBBAAAANNNNDDDDVVVV((((3333FFFF))))                                                          ____BBBBAAAANNNNDDDDVVVV((((3333FFFF))))



        2.0E0 if the eigenvalues are in descending order.  If the subroutine
     is being used to solve systems of linear equations, E21 should be set to
     1.0E0 if the coefficient matrix is symmetric and to -1.0E0 if not.

     MMMM is the number of specified eigenvalues or the number of systems of
     linear equations.

     WWWW contains the M eigenvalues in ascending or descending order. If the
     subroutine is being used to solve systems of linear equations (A-
     W(R)*I)*X(R)=B(R), where I is the identity matrix, W(R) should be set
     accordingly, for R=1,2,...,M.

     ZZZZ contains the constant matrix columns (B(R),R=1,2,...,M), if the
     subroutine is used to solve systems of linear equations.

     NNNNVVVV must be set to the dimension of the array parameter RV as declared in
     the calling program dimension statement.  On OUTPUT

     AAAA and W are unaltered.

     ZZZZ contains the associated set of orthogonal eigenvectors. Any vector
     which fails to converge is set to zero.  If the subroutine is used to
     solve systems of linear equations, Z contains the solution matrix columns
     (X(R),R=1,2,...,M).

     IIIIEEEERRRRRRRR is set to Zero       for normal return, -R         if the
     eigenvector corresponding to the R-th
        eigenvalue fails to converge, or if the R-th
        system of linear equations is nearly singular.

     RRRRVVVV and RV6 are temporary storage arrays.  Note that RV is of dimension at
     least N*(2*MB-1).  If the subroutine is being used to solve systems of
     linear equations, the determinant (up to sign) of A-W(M)*I is available,
     upon return, as the product of the first N elements of RV.  Questions and
     comments should be directed to B. S. Garbow, Applied Mathematics
     Division, ARGONNE NATIONAL LABORATORY



















                                                                        PPPPaaaaggggeeee 2222